Suppression of magnetotransport in strongly disordered graphene 
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A tight-binding model with randomly fluctuating atomic positions is studied to discuss the effect 
of strong disorder in graphene. We employ a strong-disorder expansion for the transport quantities 
and find a diffusive behavior, where the conductivity is decreasing with increasing disorder. For 
sufficiently strong disorder the magnetic field drops out of the diffusion coefficient and the conduc- 
tivity. This signals a strong suppression of magnetotransport effects, a result which is consistent 
with recent experimental observations by Morozov et al. 
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PACS numbers: 81.05.Uw, 71.55.Ak, 72.10.Bg, 73.20.Jc 

Introduction: Experimental studies on graphene re- 
vealed remarkable transport properties of this material 
0, 0, There is a robust minimal conductivity in 
the absence of an external magnetic field and a quan- 
tum Hall effect in the presence of a strong magnetic field 
P, 0] • Morozov et al. found that also for weak magnetic 
fields unusual transport properties exist: The peak of 
the magnetoresistance at zero field, the hallmark of weak 
localization in two-dimensional systems jj, is strongly 
suppressed in graphene 0, H, 0]. This suppression was 
attributed to disorder caused by ripples in the graphene 
sheet. When these ripples are removed, the quasiparticle 
mobility is significantly increased and a normal magne- 
toresistance peak appears The peak was also ob- 
served in multilayers of graphene, where ripples are less 
developed because of the higher rigidity of the material. 

In this letter we study the transport properties of 
graphene in terms of a tight-binding model for quasi- 
particles on a honeycomb lattice in the presence of a ho- 
mogeneous magnetic field. Strong disorder is introduced 
by randomly fluctuating hopping rates. The study is re- 
stricted to a system without a gate (i.e. the system is at 
the Dirac point) , where the model has a chiral symmetry. 

The usual approximation by Dirac fermions is not ap- 
plicable to strong randomness. Nevertheless, we find a 
qualitative agreement of the transport behavior between 
the tight-binding model and Dirac fermions with random 
vector potential Q or random mass 0. This indicates 
that the sublattice structure, and not the chiral symme- 
try, is responsible for the suppression of magnetoresis- 
tance peak. 

Tight-binding model on a honeycomb lattice: Hopping 
of quasiparticles on a honeycomb lattice is defined by the 
Hamiltonian 



n 



E 

R,Ri 



hR,Bf<y R c R > + h.c. 



(1) 



with quasiparticle creation (annihilation) operators c R 
(cr). The hopping rate h R R i between lattice sites R and 
B! may fluctuate from bond to bond around the average 
value t due to ripples in the graphene sheet. 



A non-orthogonal basis {01,02} (cf. Fig. 1) is used 
to express the hopping rates in terms of the sublattice 
representation by writing R = r (= r + b\) if R on A (B). 
Now r can be expressed in the basis {ai, 02}, and h R R i 
in Eq. |T]) is replaced by the matrix 



H = 







d + A 



= t 



=1,2,3 



for ci 



2.:5 



(2) 

= 0, c 2 = a 2 , and c 3 = a\ + a 2 . A r . r > = 
t' r , r 5 r ' tr - Cj are random fluctuations of the hop- 
ping term with (t' r , r ) = 0. We set t = 1 and keep in mind 
that t' and all other energies of the model are measured 
now in units of the hopping rate t. 

Although curvature can have significant effects @, 
11], it will be ignored here by assuming that the lattice 
is flat on length scales that are relevant for transport 
properties. Therefore, only fluctuations of the hopping 
rates between neighboring sites on the honeycomb lattice 
are considered. 

The quasiparticle Hamiltonian H consists of a 
translational-invariant part Hq = (H) and a term that 
describes the random fluctuations. For the latter we as- 
sume that on large scales it can be approximated by on- 
site terms, provided there is only short-range disorder. 
Such a choice leads to a random vector potential, pre- 
viously proposed for ripples in graphene by Morozov et 
al. [5[ and derived for changes in the hopping due to hy- 
bridization between different orbitals in curved graphene 
sheets and due to elastic strain [6j: 



H = Hn 



■Vl<Tl + V 2 0- 2 



(3) 



with Pauli matrices <jj. We introduce an uncorrelated 
Gaussian distribution for v with mean zero and variance 
g and write v\ = vcosa, v 2 = usina with a tunable 
parameter a. Similar models, based on Dirac fermions 
with a random vector potential, have been studied in the 
literature [1, 12, [ijj]. The seminal work by Ludwig et 



al. states that a random vector potential does not renor- 
malize the variance g, and that (using a bosonization ap- 
proach in replica space) the conductivity is not affected 
by this type of disorder Q. The density of states, on 
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FIG. 1: Lattice vectors on the honeycomb lattice. The lattice 
is divided into sublattice A and B, where nearest neighbors 
of A belong to sublattice B. 



the other hand, has a power law near the Dirac point 
with exponent (tt — q)/{tt + g). This indicates that for 
strong disorder (g > tt) the density of states diverges at 
the Dirac point, whereas it vanishes for weak disorder 
(g < tt). Physical implications of this behavior can be 
discussed in terms of the Einstein relation, where the con- 
ductivity a is given as the product of density of states p 
and diffusion coefficient D as a oc pD. Then the nonzero 
conductivity at the Dirac point (minimal conductivity) 
would imply an infinite D for g < tt and D = for g > tt. 
In other words, these findings describe a transition at 
g = tt from ballistic transport directly to localization, 
without an intermediate diffusive regime. 

For small g the exponent of the density of states can 
also be evaluated in perturbation theory as l — g/n. How- 
ever, in a recent work on the density of states a disorder 
created energy scale exp(—Tr/g)/g was detected, below 



which perturbation theory breaks down [14l ]. In other 
words, the assumption of a power law is invalid at the 
Dirac point and in its vicinity. The energy scale is simi- 
lar to the Kondo scale of the Kondo effect [IH , and was 
also found independently for the conductivity by Auslen- 
der and Katsnelson [16j. It implies a nonzero density 
of states at the Dirac point. A similar result was found 
for topological disorder [Io| . In the following, it shall be 
discussed that the nonzero density of states is associated 
with spontaneous breaking of a symmetry which leads to 
diffusion of quasiparticles. 

Magnetic field: A magnetic field B perpendicular to 
the graphene plane enters the Hamiltonian through a 
Peierls phase, i.e. the hopping terms in Hq are replaced 
by 



d rr i (i?) 



r+bj 



A-dbi 



(4) 



where A is the vector potential of the magnetic field B. 
A homogeneous magnetic field creates additional terms 
to H which are periodic in one spatial direction and pro- 
portional to (Ti and (72, respectively. 

Symmetries: H is invariant under the transformation 
H -> ^f^uef" 3 . This is related to the chiral symmetry 
discussed, for instance, in Refs. 17, 3- Here we are not 



so much interested in the symmetry of the Hamiltonian 
but in that of the two-particle Green's function (2PGF) 



at the Dirac point 

K(r,r';e) = (Tr 2 [G(r,r';ie)G(r',r;-ie)}) . (5) 

Tr2 is the trace with respect to the spinor structure, 
G(ie) = (H + ie)^ 1 is the one-particle Green's func- 
tion, and (...) is the average with respect to the random 
term of the Hamiltonian. The second Green's function 
in Eq. ([5]) can be replaced by — ay,G{r' , r; ie)a^, since 
it obeys the relation G(—ie) = —(7iG(ie)ai. Within a 
linear-response approach K(r,r';e) gives the transport 
properties: The DC conductivity at zero temperature is 

E El El m 0,0 



cxe 2 Vr?K(r,0;e)| e 



(G) 



For localized states the 2PGF decays exponentially on 
the localization length for which this expression gives a 
vanishing DC conductivity. 



To analyze the 2PGF in detail, we follow Ref. [25 1 
and write it, before averaging, as a Gaussian functional 
integral with two independent Gaussian fields, a boson 
(complex) field Xrk and a fermion (Grassmann) field ^> r k 
(* = 1,2): 

Gjf (r, r'; ie)G k -k{r' ', r; ie) 



Xr'j'Xrj^rk^r'k' exp(- Sq)V^T>X 



So is a quadratic form of the four-component field 

(Xrl,Xrl,*r2,*r-2) 



(7) 



S = -i^2(/) r ■ {H. + ie) r y<j) rl (e > 0) . (8) 

r,r' 

The use of the mixed field <p r avoids an extra normaliza- 
tion factor for the integral. (The replica trick, used in 
Ref. [Io| , is an alternative to avoid this extra normaliza- 
tion factor.) The extended Hamiltonian H = diag(H, H) 
of So acts in the boson and in the fermion sector sepa- 
rately. It is invariant under the transformation 



H UHU = H with U = cxp 



7/>CT 3 \ 

^3 J 



(9) 



for Grassmann variables ip and -tp, whereas ie in Eq. ([8]) 
causes a symmetry breaking. This symmetry is central 
for the transport properties, since it determines the large- 
scale properties of the 2PGF by creating a massless mode. 
The latter is a consequence of a possible spontaneous 
symmetry breaking at the Dirac point e — * [9]. This 
indicates that transport at the Dirac point is qualita- 
tively different from transport away from the Dirac point, 
where this symmetry is explicitly broken by the Fermi en- 
ergy. An important question in this context is how the 
symmetry affects average quantities, and whether or not 



3 



spontaneous symmetry breaking occurs. The signature 
of the latter is a nonzero density of states at the Dirac 
point 26]. (The method of Ref. [26j can be applied to 
the present model to prove that the average density of 
states is nonzero.) 

Averaging the 2PGF over the Gaussian distribution 
of v r means replacing exp(— So) by (exp(— So)) on the 
right-hand side of Eq. ([7]). The average quantity can be 
written again as an exponential function (exp(— So)) = 
exp(— Si), where the new function Si contains also quar- 
tic terms of the field </>, a consequence of the Gaussian 
distribution. Then it is convenient to transform the in- 
tegration variables as 

( XrXr Xr^r\^ r . _ ( Qr ®r\ nrn 
[VrXr Mr ^6 r -iPr) ' 

where Q r , P r are symmetric 2x2 matrices, and r , 
G> are 2x2 matrices whose elements are independent 
Grassmann variables. The average 2PGF now reads as a 
correlation function in the new field Q r : 

K(r,r';e) = ~ [ Tr 2 (e r , 7 )Tr 2 (e r 7) exp(-S 2 )P[Q] 



S 2 = V -Trg(Q 2 ) - ln[detg[ff + ie - 2 7 Q]] . (12) 
,9 

Trg is the graded trace, detg the graded determinant [25| . 
and 7 is a 4 x 4 matrix, consisting of the linear combi- 
nation 71 cos a + 72 sin a, where 7/ = diag(aj, <7j). The 
integration over Q can be performed in a saddle-point 
(SP) approximation, where the SP satisfies the extremal 
condition <5qS 2 = 0, which reads explicitly 

Q s r p = 2g[H + 2 7 Q SP ]- 1 7 . (13) 

A consequence of the symmetry ([9]) is that with Q SP also 
is a solution of the SP equation at the Dirac 
point e = 0. 

The right-hand side of Eq. (fTB"]) is translational invari- 
ant in the absence as well as in the presence of a homo- 
geneous magnetic field. This is a consequence of the fact 
that the diagonal Green's function G rr describes closed 
loops of quasiparticles. These loops depend only on the 
flux which they enclose. The translational-invariant G rr 
has a constant SP solution Q SP = —1(77/2)7 which sat- 
isfies 

r) = Aig(H a + ie + iry)" 1 . (14) 

X] can be interpreted as a self energy of the average one- 
particle Green's function (Ho + ie + irj)^ 1 20]. Then the 
average density of states at the Dirac point is p = rj/g. 
A SP solution 77 > exists even in the limit e — * 0, 
which reflects spontaneous symmetry breaking. rj in- 
creases with g and reaches for large g the asymptotic 



regime with rj ~ 2y / g, where the density of states decays 
like p ~ .g -1 / 2 . 

The invariance of the SP equation requires the integra- 
tion over all SPs and their vicinities. This leads to an SP 
manifold, generated by the symmetry transformation in 
form of the nonlinear field 

Q' r = -i|u+ 7 U r = -i| 7 U? . (15) 

U r is the matrix U of Eq. ([9]) with space-dependent 
Grassmann fields ip r and ip r . The integration with re- 
spect to the chiral symmetry is not taken into account 
here because the corresponding fields are perpendicular 
to the Grassmann fields in leading order. Thus they do 
not contribute to the average 2PGF. 

Replacing Q r by QJ. in S 2 provides two major simpli- 
fications: (1) the first term in expression (fT2"|) vanishes, 
since QJ 2 is proportional to the 4x4 unit matrix. (2) 
the second term becomes 

-lndctg(U- 1 (Ho + ie)U- 1 + i 7 y ), (16) 

since detgU = 1. This expression can be expanded ei- 
ther in powers of 77 -1 or in powers of 77. The problem of 
the latter case is that at low energies, which are relevant 
at the Dirac point, the expansion terms (Ho + ie)~ l are 
arbitrarily large. Consequently, this expansion cannot 
be controlled. The expansion in powers of ?7 _1 , on the 
other hand, has always small terms for a smoothly vary- 
ing U" 1 , provided r] is not too small in comparison with 
energies of these modes. The expansion in powers of r? -1 
can also be justified by using large values of 77 formally 
and then continue it down to physically reasonable val- 
ues of 77. Since no singularity appears, this extrapolation 
is at least qualitatively correct. Thus the second-order 
term is a good approximation for low energies. It gives 
for the expression in Eq. (|16p 

S 2 = V(e^r' + CV r ,)VvVv + o(tT 3 ) , (17) 

77 * — ' 

r,r' 

where C rr i is 

— ( Tr 2 [Ho t rr'Ho,r'f] $r,r' — Tr 2 [Ho,rr' Ho.r'r] 1 ■ 
' f 

(18) 

The average 2PGF of Eq. (IT]) then reads 

K(r,r';e)^^(e + C)-] . (19) 

This is the main result of our SP calculation. It means 
that details of the transport properties depend only on 
the SP solution 77 and the average Hamiltonian Ho in the 
correlation C. A similar expression was found for the 
2PGF of weakly disordered Dirac fermions [9j. 
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Discussion: 1) B = 0: Translational invariance of C r 
suggests to study its Fourier components 



C(q) 



cos qi — cos 



'{V3 qi +q 2 )/2 



(20) 



or for the more general case (e.g. for a low-energy ap- 
proximation of Hq by Dirac fermions), the Fourier com- 
ponents 

C ^ = 2~7? / [ Tl " 2 ^^ ~~ Tl " 2 ( h k+q/2 h k-q/2)] ■ 

Here /i^ are the Fourier components of H , i.e. for Dirac 
fermions we have hk = kio~i + fc2<72- C(<?) vanishes at 
g = and describes the dispersion of the two-particle 
excitations. In constrast to quasiparticles with Hamil- 
tonian Hq, which have low-energy excitations near the 
Dirac nodes k ^ 0, the two-particle excitations at low 
energies appear around q = 0. This implies that K{q; e) 
is a diffusion propagator with diffusion tensor 



_ 1 d 2 C(q) 
13 2 dqidqj 



= 2 f Tr Jdhkdhk 



(21) 



Then the main contribution to diffusion comes from the 
vicinities of the Dirac points, where the slope of hk is 
the strongest. This would justify an approximation of 
hk by Dirac fermions. Furthermore, in agreement with 
physical intuition, D is monotoneously decreasing with 
increasing disorder. Finally, the DC conductivity a can 
be evaluated via the Kubo formula as 



a 2 



(22) 



in units of e 2 /h. This result indicates that the minimal 
conductivity is not constant but decreases with increasing 
disorder. For the regime of strong disorder the conduc- 
tivity is proportional to <? _1 , which was also found for 
topological disorder (loj . 

2) B 0: The Hamiltonian Ho appears in C rr > only 

as 

Tr 2 [H 0trr ,H 0ir , r ] = \d rr ,{B)\ 2 + \d r , r {B)\ 2 . 
For a given pair of nearest neighbors (r, r f ) we have only 



i(j> r 



>8 / A. I 2 
u r f ,r±Cj | 5 



one contribution for each pair, namely |e 
such that the Peierls phase drops out: \d rr '(B)\ 2 — 
|d rr /(0)| 2 . Thus, the magnetic field dependence of the 
average 2PGF in Eq. P^|) is due to the SP solution r/. 
The latter can be evaluated from Eq. ([14")) by expanding 
the right-hand side in powers of ry _1 : 

>) 2 = 4j-iE Tr 2 [Ho.rr'Hoyr] + o{n- 3 ) . (23) 

r' 

Again, the B dependence appears only in higher orders 
of the expansion. While rj increases monotoneously with 
g, its B dependence is suppressed. 



Disorder due to potential scattering by impurities has 
been ignored in our approach. Including it in the calcu- 
lation would affect transport properties at low tempera- 
tures significantly, since it causes localization (2t| ■ It is 
not clear, however, whether or not this type of disorder 
is relevant in graphene at the Dirac point. 

We conclude that quasiparticles on the honeycomb lat- 
tice with random bonds are subject to diffusion. The 
diffusion coefficient as well as the conductivity always 
decrease with increasing disorder. For sufficiently strong 
disorder the magnetic-field dependence drops out of the 
transport quantities. This reflects the suppression of an 
external magnetic field by disorder, which is consistent 
with experimental observations of magnetotransport in 
graphene [1,0]. 

I would like to thank Balazs Dora, Andre Geim, 
Mikhael Katsnelson, and Maria Vozmediano for stimu- 
lating discussions. 
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